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Abstract 

We demonstrate theoretically and numerically that the fundamen¬ 
tal Raman soliton of the widely used nonlinear Schroedinger equation 
(NLSE) with a linear approximation of the Raman gain [reduced NLSE) 
is metastable. It can propagate for hundreds of dispersion lengths along 
the optical fibre before eventually disappearing due to a peculiar insta¬ 
bility, leading to a collapse. The noise eigenfunction analysis agrees well 
with the results obtained via direct pulse propagation simulations. This 
instability is not present when modelling the Raman effect via a full con¬ 
volution, and thus the reduced NLSE often leads to unphysical results, 
and should be avoided. 

The dynamics of a short and intense pulse in optical fibres is determined by 
several physical effects. One major effect is the stimulated Raman scattering 
(SRS) which becomes particularly important for sub-picosecond pulses [HIS]- 
In the presence of the Raman effect combined with Kerr nonlinearity and in 
the absence of higher-order dispersion (HOD), Raman states with a single or 
double peak structure were numerically found by solving the reduced nonlinear 
Schroedinger equation (NLSE), where a linear approximation of the Raman 
gain is used [see Eqs. Q - ([^ below], by Akhmediev et a/.[3]- In Ref. [5], the 
Raman two-peak states were also shown to be metastable in the sense that they 
can propagate for long distances before eventually disappearing due to internal 
collapse. A few years ago, in a series of experiments gna, short and intense 
pulses have been launched in highly nonlinear photonic crystal fibres (PCFs, 
see also 0 )- The formation of long-lived, two-peak soliton states for specihc 
’magic’ input pulse energies was observed. In these experiments, it was clear 
that this phenomenon is quite universal for a large variety of highly nonlinear 
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optical fibres, especially in the deep anomalous group-velocity dispersion (GVD) 
regime far from zero-GVD points, where the impact of HOD is minimised. Such 
localised states always show peaks of differing amplitude with a well-defined 
peak power ratio r. Moreover, the temporal separation between the two peaks 
is also uniquely determined by r and the Raman parameter A full theoretical 
understanding of the excitation of Raman multipeak states (with an arbitrary 
number of peaks) in solid-core PGFs was provided in Ref. [71 [H]. Based on 
a ’gravity-like’ potential approach, in Ref. |7] we derived simple equations for 
the ’magic’ peak power ratio r and the temporal separation between pulses 
forming these multipeak states. We also developed a model to calculate the 
magic input power of the input pulse around which the phenomenon can be 
observed. In Ref. [71 [S] we also numerically demonstrated that all Raman 
multipeak states are metastable, and they become more and more unstable when 
increasing the number of peaks. This instability occurs in both the reduced and 
the full NLSE, where in the latter case the Raman response function is used 
without approximations. 

The physical metastability of Raman multipeak states has been thus proved 
unambiguously in our previous work |7]. However, the stability of the single 
Raman peak state (further referred to as the fundamental Raman soliton, FRS 
in brief) has not yet been determined conclusively. In Ref. |3], the FRS was 
claimed to be stable via numerical simulations of the reduced NLSE for a prop¬ 
agation distance up to about 100 dispersion lengths. However, as we will show 
later, this propagation distance is not sufficient to prove that the FRS is stable 
in the reduced model. In Ref. [3], the authors used arguments based on the 
Evans function, concluding again that the FRS of the reduced NLSE is stable. 
However, we show here that this is not the case. 

Our analysis explains and confirms previous numerical simulations performed 
by Erkintalo et al. nni and Solli et al. E] , who observed pulse collapse and energy 
non-conservation when using the reduced NLSE. In this latter paper, famous for 
the first experimental observation of rogue waves in optics, the authors remark 
that “The rogue waves have a number of other intriguing properties warranting 
further study. For example, they propagate without noticeable broadening for 
some time, but have a finite, seemingly unpredictable lifetime before they sud¬ 
denly collapse owing to cumulative effects of Raman scattering. [...] The decay 
parallels the unpredictable lifetimes of oceanic rogue waves”. It is shown here 
that this collapse is totally unphysical and is a byproduct of the instabilities 
introduced by the linearized Raman gain, as hinted in m ■ 

In this Letter we theoretically and numerically demonstrate that the FRS 
of the reduced NLSE is unstable: it can propagate for several hundreds of dis¬ 
persion lengths along the fibre without signihcant changes in its profile, but 
eventually it will collapse due to a peculiar instability. However, long-distance 
simulations of the full NLSE conclusively show that the Raman soliton is per¬ 
fectly stable in this physically more precise model. This underlines important 
physical inconsistencies of the reduced NLSE, which should therefore be care¬ 
fully avoided for numerical modelling. 

For wide enough pulses (temporal width > 100 fs) the generalised NLSE 
governing the pulse propagation in silica optical fibres and solid-core PGFs can 
be simplihed in dimensionless units in the laboratory frame {z,t) |3], by using a 
linear approximation of the Raman response function, giving rise to the reduced 
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NLSE: 


( 1 ) 


'idzip + = 0 , 

where t) is the electric field envelope rescaled with the fundamental soliton 
power Ps = 1/3217(7^0)) 1^2 is the GVD coefficient at the reference frequency, and 
7 is the fibre nonlinear coefficient. Variables 2 and t are dimensionless space and 
time, respectively, rescaled with the second order dispersion length Pd = to/l/^2| 
and with the input pulse duration to- The last term of Eq. Q, responsible for 
the Raman effect, produces a constant Raman soliton self-frequency shift (SSFS) 
in silica fibres [laiin!, Tji = Tji/to being a small parameter, where Tr is the 
Raman response time, approximately equal to 2 fs in silica. Solitons subject to 
the Raman effect shift continuously towards the red part of the spectrum due 
to SSFS, leading to a constant (negative) acceleration of the soliton in the time 
domain [iiin]. 

In the reference frame of an intense accelerating soliton, one can operate the 
so-called Gagnon-Belanger gauge transformation t) = f{x, exp \iz {q — b'^z 
where x = z, ^ = t — hz^ 12 , and b = 32TRg^/15. This leads to the following 
equation for f{x,^) [CTI^: 

idxf + ^5|/ - {q + b^)f + \f\^f - TR/^jd/p) = 0, (2) 

where q is the wavenumber of the most intense soliton, proportional to its peak 
power. Note that in mm the equivalent form of Eq. ([^ does not contain 
the first term idxf because in these works only the profile of the stationary 
solution was considered. To this aim, one can safely ignore the space derivative 
dxf to obtain an ODE as done in HE]. However, in this work, we also want to 
investigate the dynamics of the FRS in the accelerated frame (x, ^), and thus we 
wish to keep the space derivative dxf in Eq. ([^. If one neglects all the nonlinear 
terms, Eq. ([^ would correspond to the classical Schroedinger equation for a 
unitary mass particle of energy q subject to a ’gravitational’ potential U{^) = b^ 

[TSlfTB] . 

Equation ([^ is non-integrable, and its localised solutions cannot be ana¬ 
lytically obtained. Moreover, the soliton is not stationary and the traditional 
version of the Vakhitov-Kolokolov criterion m and the Evans function approach 
H seem not to lead to the right conclusions on the stability. 

The profile of the stationary localised solutions to Eq. (|^ was numerically 
obtained in HE]- These localised solutions, theoretically, do not change when 
the space variable x evolves, i.e., dxf = 0. However, if these solutions are not 
stable, their profiles can still change significantly during propagation as shown 
later. An example of the FRS profile as a function of the delay ^ is shown in Fig. 
[^a) (red dashed-dotted curve labeled fs), for a soliton wavenumber q = 0.1. 
This FRS shown in Fig. Ba) is a real function, and in what follows we will always 
assume that fs is real without losing generality. This is due to the properties of 
Eq. ([^ in which if fs is a solution then fsexp(ia) is also a solution with a being 
an arbitrary constant phase. For any value of the soliton wavenumber q, such 
solutions always show a characteristic Airy tail with oscillations on the leading 
edge of the pulse due to tunnelling in the linearized Schroedinger equation. The 
Airy tail is quite small when using the physically relevant parameters, and thus, 
in Fig. j^a) one cannot see it without strongly zooming the leading edge. Airy 
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tails will be more pronounced if the slope b of the gravity-like potential gets 
larger, due to the exponential increase of the tunneling probability. 

Now we use the linear stability approach to analyse the stability of the FRS. 
To this end we insert the Ansatz f = fs+g into Eq. ([^, where the complex field 
g is a weak noise added to the FRS /g. After linearising Eq. with respect 
to g and splitting g into real and imaginary parts {g = gi + * 32 ) one gets the 
following system of linear equations for the two real fields gi and ( 72 : 

- d^g 2 + ^d^gi -{q + b^gi + 3gif^ - Tagid^if^) - 2Tiif^d^{f^gi) = 0, (3) 

dxgi + - (9 + b^)g2 + g2fs - TRg2d^if^) = 0. (4) 

It is clear from Eqs. © and 0 that gi and (72 are coupled during propaga¬ 
tion. However, we first study the stationary solutions of gi^ 2 -> so one can ignore 
the terms dxg 2 and dxgi in Eqs. ^ and Q. In this case, one obtains two 
following independent linear equations for the noise eigenfunctions: 


^{51 + Pigi + ‘2Tiifsd^gi — 0 , 

(5) 

— ^d\g2 + ^252 = 0 , 

( 6 ) 


where the potentials felt by gi ^2 are respectively Pi = q + b^ — 3/f -I- 
and P 2 = q + b^ — + 2TRfsd^fs. Note that Eqs. 0 and 0 are formally 

identical to two stationary Schroedinger equations with potentials Pi ^ 2 , but the 
third term in Eq. 0 provides gain for gi , leading to a growth for the real part 
of the noise - and thus also to the imaginary part (72 via the dynamical coupling 
along the propagation. This will lead to instabilities of the FRS, and is the 
main result of this Letter. 

In Fig. 0a) we plot Pi ^2 (black dashed and orange solid curves, respectively) 
as functions of the delay ^ for the specific FRS /g. 




Figure 1: (Color online) (a) The profile of the fundamental Raman soliton fs with q 
= 0.1 (red dashed-dotted curve) and the two potentials Pi, 2 (black dashed and orange 
solid curves, respectively) as functions of the delay (b) The profiles of the two noise 
eigenfunctions gi (green dashed curve) and 772 (red solid curve) corresponding to the 
two potentials Pi ,2 in (a), respectively. < 71,2 are normalised such that the maximum 
values of \gi, 2 \ are equal to the maximum value of the fundamental Raman soliton /a. 
Note that the normalised <72 is identical to /a. Raman parameter is tr = 0.1. 
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Now we find the profiles of the two noise eigenfunctions 31,2- By numerically 
solving Eqs. ([^ and Q with the given fg, we are able to obtain the localised 
solutions for gi 2 as shown in Fig. Bb) (green dashed curve and red solid 
curve, respectively). For comparison we also plot the FRS fs in Fig. Sb) (the 
blue curve hidden behind the red solid curve). As mentioned above, the noise 
eigenfunctions gi^2 must be much weaker than the FRS fs- Just for the sake 
of comparison, the eigenfunctions 51^2 plotted in Fig. Sb) are normalised such 
that the maximum values of |5i,2| are equal to the maximum value of the FRS 
fs- In Fig. [^b) one can easily see that the functional shape of (72 is identical to 
fs, therefore, the red and blue curves representing these two functions are on top 
of each other in Fig. Sb). This is not surprising because if one replaces g2 in 
Ep. (joj) by fs^ then one gets immediately the stationary hip. ([|) (for id^f = 0 ). 

In what follows we will compare the noise eigenfunctions gi^2 (which are 
theoretically obtained above) with the noise generated during propagation of 
the FRS. Figure l^a) shows the propagation of the FRS [whose initial profile is 
fs as plotted in Fig. [^a)] in the accelerated reference frame (a;,^). As clearly 
seen in Fig. [^a), when the propagation distance x < 1000 the profile of the 
FRS does not change significantly during propagation. However, for larger 
propagation distances, the noise grows noticeably, then at a distance x ~ 1500 
the FRS is completely destroyed by a powerful instability and one only gets 
an approximately uniform noisy background. We want to emphasise that this 
behaviour of the FRS during propagation is not an artefact due to the numerical 
inaccuracy accumulated during the propagation, because one always gets this 
instability even when the numerical accuracy is greatly enhanced (larger delay 
window [fniinjimax], Smaller discretisation step for both the and a;-axis). So, 
Fig. I^a) is a good evidence that, in the reduced NFS, the FRS is metastable like 
all other Raman multipeak states - it can propagate for long distances without 
significantly changing its profile, but finally it will collapse for a specific critical 
distance. 

Now we define the deviation of the FRS profile at the current value of the 
propagation distance x from the initial FRS fs profile during its propagation as 
follows: 

d{x,f) = di{x,f) + id2{x,f) = f{x,f,) - fs{ 0 , <J) 

where f{x, ^) is the profile of the FRS at the distance x that one gets by numer¬ 
ically solving Eq. ([^, with the input field being fs, whereas di and d2 are the 
real and imaginary parts of d, respectively. These deviation components di_2 can 
be interpreted as the real and imaginary parts of the noise during propagation 
of the FRS. In a stable situation, as mentioned above, this deviation must be 
equal to zero, i.e., d{x,f,) = di{x,f^) = d2{x,ff) = 0 . However, since the FRS is 
not stable, this deviation will differ from zero and grow during propagation. In 
Figs. ^h,c ,d) we plot the real (di - green curve) and imaginary (c?2 - blue curve) 
parts of d{x, for three values of the propagation distance, x = 5 , 50 , and 100 , 
respectively. The profiles of the noise eigenfunctions gi^2 are also plotted in 
Figs. [^b,c,d) with black and red curves, respectively. 51,2 in Figs. [^b,c,d) are 
re-normalised such that max(gii) = max(di), and min((72) = min(d2). At the 
beginning of the propagation when x = 5 [Fig. [^b)], the theoretically obtained 
eigenfunctions 2 differ significantly from the corresponding deviation compo¬ 
nents di^2- However, for larger propagation distances, as clearly shown in Fig. 
[|;c) {x = 50 ) and Fig. [|;d) (x = 100), the deviation components of the noise 
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Figure 2: (Color online) (a) Propagation of the fundamental Raman soliton and its 
instability. (b,c,d) The profiles of the real (di - green curve) and imaginary parts of 
the deviation (d 2 - blue curve) of the FRS profile during propagation as compared to 
the FRS profile /s at the input for three values of the propagation distance a; = 5, 50, 
and 100, respectively. The profiles of re-normalised eigenfunctions g \^2 are also plotted 
in (b,c,d). Fig. 2(a) is on a logarithmic scale 


di ^2 have profiles which are almost identical to gi, 2 , fully confirming our theo¬ 
retical analysis based on Eqs. (3]|4). Note that when the propagation distance is 
increased, the noise gradually grows in strength, for instance, when a: = 5 [Fig. 
I^b)] max(di) ~ 2.10“®, whereas when x = 100 [Fig. [^d)] max(di) ~ 3.10“"^. 
The profile difference between di 2 and 51^2 in Fig. Sb) is understandable, be¬ 
cause the deviation components di ^2 of the noise during propagation need time 
to build up and acquire the profiles of their corresponding eigenfunctions gi^ 2 - 
For much larger propagation distances {x > 1000), as shown in Fig. [§&), the 
noise grows significantly, and as a result the condition that the noise must be 
much weaker than the FRS is not longer satisfied. Evidently, in this case, the 
procedure for obtaining the gi ^2 discussed above is not valid. 

In Fig. I^a) we show the propagation of a FRS in the laboratory frame 
(z,t) when using the full Raman convolution, see for instance Eq. (2.3.36) of 
Ref. [T]. We use 3000 frames along the z-axis to show results in Eig. §a). 
After reaching the bottom part of the time window, the soliton will re-emerge 
in the upper part due to the periodic boundary conditions. However, since 
there is never self-interaction and the time windows is large, this simulates 
an undisturbed and unlimited propagation. We have propagated the FRS for 
6000 dispersion lengths and concluded that there is no instability that leads to 
the soliton collapse. In Fig. [^c) we show the spectral evolution of the FRS 
during propagation. The temporal and spectral profiles of the FRS are plotted 
in Figs. ib) and (d), respectively, for different propagation distances. As 
clearly shown in Fig. [^ both temporal and spectral profiles of the FRS are well 
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Figure 3: (Color online) (a,c) Temporal and spectral evolution of a FRS using the 
full Raman convolution, respectively. (b,d) The temporal and spectral profiles of the 
FRS, respectively, for different propagation distances 2 . 


conserved during propagation. It is important to note that the input condition 
for simulation results in Fig. [^a) is the numerically obtained fs of the reduced 
NLSE shown in Fig. [^a), which is not the exact soliton solution of the NLSE 
with the full Raman convolution. One can thus conclude that the physical FRS 
is perfectly stable. 

In conclusion, we have investigated the stability issue of the fundamental 
Raman soliton of the reduced nonlinear Schroedinger equation in optical fibres, 
where the Raman gain is approximated as linear. We demonstrate theoretically 
that the FRS is metastable in the reduced NLSE, in the sense that it can 
propagate several hundreds, even thousands of dispersion lengths along the fibre 
without significant changes in its prohle. However, with larger propagation 
distances the noise is amplified, and finally, the structure of the fundamental 
Raman soliton completely collapses due to a powerful instability, confirming and 
explaining previous numerical observations mm- The eigenfunctions of the 
noise obtained via our simple theoretical approach agree very well with the noise 
of the FRS obtained directly through pulse propagation simulations. However, 
we show that the FRS of the full NLSE is perfectly stable. This underlines 
the dangers of using approximations of the Raman response functions in any 
dynamical simulations. 

This work is supported by the German Max Planck Society for the Advance¬ 
ment of Science (MPG). F.B. would like to thank M. Erkintalo, J. Dudley and 
G. Genty for useful comments. 
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